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Abstract 

We investigate the approach to stable and metastable equilibrium in Ising models using a cluster 
representation. The distribution of nucleation times is determined using the Metropolis algorithm 
and the corresponding (j) 4 model using Langevin dynamics. We find that the nucleation rate is 
suppressed at early times even after global variables such as the magnetization and energy have 
apparently reached their time independent values. The mean number of clusters whose size is 
comparable to the size of the nucleating droplet becomes time independent at about the same 
time that the nucleation rate reaches its constant value. We also find subtle structural differences 
between the nucleating droplets formed before and after apparent metastable equilibrium has been 
established. 
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I. INTRODUCTION 



Understanding nucleation is important in fields as diverse as materials science, biological 
physics, and meteorology . Fundamental progress was made when 

Gibbs assumed that the nucleating droplet can be considered to be a fluctuation about 
metastable equilibrium, and hence the probability of a nucleating droplet is independent of 



time 



lOj. Langer ll| has shown that the probability of a nucleating droplet can be related 



to the analytic continuation of the stable state free energy in the limit that the metastable 
state lifetime approaches infinity. Hence the assumption by Gibbs is valid in this limit. It has 
also been shown that the Gibbs assumption is correct in systems for which the interaction 



range R — > oo 



13|- 



For metastable states with finite lifetimes equilibrium is never reached because a large 
enough fluctuation would initiate the transformation to the stable state. However, if the 
probability of such a fluctuation is sufficiently small, it is possible that systems investigated 
by simulations and experiments can be well approximated as being in equilibrium. Hence, 
for metastable lifetimes that are very long, we expect the Gibbs assumption to be a good 
approximation. 

In practice, nucleation is not usually observed when the lifetime of the metastable state 
is very long. Processes such as alloy formation, decay of the false vacuum, and protein 
crystallization generally occur during a continuous quench of a control parameter such as 
the temperature. It is natural to ask if the nucleation process that is observed occurs when 
the system can be reasonably approximated by one in metastable equilibrium. If so, the 
nucleation rate will be independent of time. 

It is usually assumed that metastable equilibrium is a good approximation when the 
mean value of the order parameter and various global quantities are no longer changing 
with time. As an example, we consider the nearest-neighbor Ising model on a square lattice 
and equilibrate the system at temperature T = 4T c /9 in a magnetic field h = 0.44. The 
relatively small value of the linear dimension L = 200 was chosen in order to avoid nucleation 
occurring too quickly. At time t — the sign of the magnetic field is reversed. In Fig. [T] we 
plot the evolution of the magnetization m(t) and the energy e(t) per spin using the Metropolis 
algorithm. The solid lines are the fits to an exponential function with the relaxation time 
T g m 1.5. In the following we will measure the time in terms of Monte Carlo steps per spin. 
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(a) m(t). (b) e(t). 

FIG. 1: The evolution of the magnetization m(t) and the energy e(i) per spin of the nearest- 
neighbor Ising model on a square lattice with linear dimension L = 200 using the Metropolis 
algorithm. The system was prepared at temperature T = 4T c /9 in the external magnetic field 
h = 0.44. At time t = the sign of the magnetic field is reversed. The solid lines are fits to 
an exponential function with relaxation time r g = 1.5 and 1.2 respectively. (Time is measured in 
Monte Carlo steps per spin.) The data is averaged over 5000 runs. 



A major goal of our work is to address the question, "Can the system be treated as being 
in metastable equilibrium for t ^ r 9 ?" 

If the nucleation rate is independent of time, the probability of a nucleating droplet 
occurring at time t after the change of magnetic field is an exponentially decreasing function 
of time. To understand this dependence we divide the time into intervals At and write the 
probability that the system nucleates in a time interval At as XAt, where the nucleation 
rate A is a constant. The probability that nucleation occurs in the time interval (N + 1) is 
given by 

P N = (1 - XAtfXAt. (1) 
If we assume that A At is small and write N = tj At, we can write 

P(t)At = (1 - AAt) </A *AAt -» e" A 'AAt, (2) 

where P(t)At is the probability that the system nucleates at a time between t and t+ At after 
the change of the magnetic field. In the following we ask if the nucleation rate and the mean 
values of the order parameter and other thermodynamic quantities become independent 
of time at approximately the same time after a quench or is the approach to metastable 
equilibrium more complicated? 
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In Sec. |TT] we determine the probability distribution of the nucleation times and find that 
the nucleation rate becomes a constant only after a time T nequ ii that is much longer than the 
relaxation time r g of m(t) and e(t). In Sec. II III we study the microscopic behavior of the 
system and determine the relaxation time r s for n s , the mean number of clusters of size s, 
to approach its equilibrium value [jjjl . Our main result is that t s is an increasing function of 
s, and the time required for n s to reach its equilibrium value is the same order of magnitude 
as T noqu ii for values of s comparable to the nucleating droplet. That is, the time for the 
number of clusters that are the size of the nucleating droplet to reach its equilibrium value 
is considerably longer than the time for the mean value of the order parameter to become 
independent of time within the accuracy that we can determine. 

In Sees. IIVI and IVl we show that there are subtle differences between the structure of the 
nucleating droplets which occur before and after metastable equilibrium appears to have 
been achieved. This difference suggests the possibility of finding even greater differences in 
the nucleating droplets in systems of physical and technological importance. We summarize 
and discuss our results in Sec. IVI1 In the Appendix we study the evolution of the clusters 
after a quench to the critical temperature of the Ising model and again find that that the 
clusters equilibrate in size order, with the smaller clusters equilibrating first. Hence in 
principle, an infinite system will never equilibrate. How close to equilibrium a system needs 
to be and on what spatial scale so that it can be treated by equilibrium methods depends 
on the physical process of interest. 

II. DISTRIBUTION OF NUCLEATION TIMES 

We simulate the Ising model on a square lattice with interaction range R with the Hamil- 
tonian 



where h is the external field. The notation <i,j> in the first sum means that the distance 
between spins % and j is within the interaction range R. We studied both nearest-neighbor 
(R = 1) and long-range interactions (R > 20). The interaction strength J is scaled as 
J = 4/q, where q = 2R(R + 1) is the number of interaction neighbors per spin. The 
external field h and the temperature are measured in terms of J. All of our simulations 
are at temperature T = 4T c /9, where T c is the critical temperature. For R = 1 the critical 




(3) 



<ij> 
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temperature is T c « 2.269. For R > 20 the mean field result T c = 4 is a good approximation 



to the exact value of the critical temperature [21|. As discussed in Sec. [I] the system is 
equilibrated in a magnetic field h. The time £ = corresponds to the time immediately after 
the magnetic field is reversed. 

The clusters in the Ising model are defined rigorously by a mapping of the Ising critical 



. Two 



point onto the percolation transition of a properly chosen percolation model 
parallel spins that are within the interaction range R are connected only if there is a bond 
between them. The bonds are assigned with the probability p\> = 1 — e~ 2/3J for R = 1 and 
Pb = 1 — e~ 2 P J<yl ~ pS} near the spinodal, where p is the density of the stable spins, and (3 is the 
inverse temperature. Spins that are connected by bonds form a cluster. 

Because the intervention method [jjj] of identifying the nucleating droplet is time con- 
suming (see Sec. II Vj) . we use a simpler criterion in this section to estimate the nucleation 
time. We monitor the size of the largest cluster (averaged over 20 bond realizations) and 
estimate the nucleation time as the time when the largest cluster first reaches a threshold 
size s* . The threshold size s* is chosen so that the largest cluster begins to grow rapidly 
once its size is greater than or equal to s*. Because s* is larger than the actual size of the nu- 
cleating droplet, the nucleation time that we estimate by this criterion will be 1 to 2 Monte 
Carlo steps per spin later than the nucleation time determined by the intervention method. 
Although the distribution function P(t) is shifted to slightly later times, the nucleation rate 
is found to be insensitive to the choice of the threshold. 

Figure [2] shows P(t) for R = 1 and h = 0.44, where P(t)At is the probability that 
nucleation has occurred between time t and t + At. The results for P(t) were averaged over 
5000 runs. The mean size of the nucleating droplet is estimated to be approximately 25 
spins for this value of h. Note that P(t) is an increasing function of t for early times, reaches 
a maximum at t — r noqui i « 60, and fits to the expected exponential form for t > T nequ Q. 
The fact that P(t) falls below the expected exponential for t < r neq uii indicates that the 
nucleation rate is reduced from its equilibrium value and that the system is not in metastable 



12] and 



equilibrium. Similar nonequilibrium effects have been observed in Ising-like 16 
continuous systems [3]. We conclude that the time for the nucleation rate to become 
independent of the time after the change of magnetic field is much longer than the relaxation 
time T g ~ 1.5 of the magnetization and energy. We will refer to nucleation that occurs before 
metastable equilibrium has been reached as transient nucleation. 
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(a) P(t). (b) lnP(i). 

FIG. 2: The distribution of nucleation times P(t) averaged over 5000 runs for the same system as 
in Fig. [TJ The threshold size was chosen to be s* = 30. (The mean size of the nucleating droplet 
is ~ 25 spins.) (a) P(t) begins to decay exponentially at r nequ ii ~ 60. The nucleation rate after 
equilibrium has been established is determined from the log-linear plot in (b) and is A ~ 9 x 10 
(see Eq. ©). 



In order to see if the same qualitative behavior holds near the pseudospinodal, we simu- 
lated the long-range Ising model with R = 20 and h = 1.258. In the mean-field limit R — > oo 
the spinodal field is at h s = 1.2704 (for T = 4T c /9). A plot of m(t) for this system is shown 



in Fig. 3(a) and is seen to have the same qualitative behavior as in Fig. [2] for R = 1; the 



relaxation time r g m 4.5. In Fig. 3(b) the distribution of nucleation times is shown, and we 



see that P(t) does not decay exponentially until t ^ r neq uii = 40. According to Ref. LL9j, r neq uii 
should become comparable to r g in the limit R — > oo because the free energy is described 
only by the magnetization in the mean- field limit. We find that the difference between r neq uii 



and T g is smaller for R = 20 than for R = 1, consistent with Ref. 
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III. RELAXATION OF CLUSTERS TO METASTABLE EQUILIBRIUM 

Given that there is a significant time delay between the relaxation of the magnetization 
and the energy and the equilibration of the system as measured by the nucleation rate, 
it is interesting to monitor the time-dependence of the cluster-size distribution after the 
reverse of the magnetic field. After the change the system gradually relaxes to metastable 
equilibrium by forming clusters of spins in the stable direction. How long is required for 
the number of clusters of size s to reach equilibrium? In particular, we are interested in the 

6 




(a) m(t). (b) HP(t)). 

FIG. 3: (a) The evolution of m{t) for the long-range Ising model on a square lattice with R = 20, 
/i = 1.258, and L = 500. The solid line is an exponential fit with the relaxation time T g ~ 4.5. 
The data is averaged over 2000 runs, (b) The distribution of nucleation times P(t) for the same 
system and number of runs. P(t) decays exponentially for t > r nequ ii ~ 40. The nucleation rate 
once equilibrium has been established is A = 6.4 x 10 -2 . The mean size of the nucleating droplet 
is rj 300 spins. 

time required for clusters that are comparable in size to the nucleating droplet. 

We first consider R = 1 and monitor the number of clusters n s of size s at time t. To 
obtain good statistics we chose L = 200 and averaged over 5000 runs. Figure H] shows the 

0.020 1 1 1 1 1 1 




FIG. 4: The evolution of the number of clusters of size s = 6 averaged over 5000 runs for R = 1 
and the same conditions as in Fig. [TJ The fit is to the exponential form in Eq. @ with r s ~ 8.1 
and n Sj00 = 0.0175. 
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FIG. 5: (a) The equilibration time t s as a function of the cluster size s for R = 1 and h = 0.44 the 
same conditions as in Fig. [TJ The s-dependence of t s is approximately linear. The extrapolated 
value of t s corresponding to the mean size of the nucleating droplet 25 spins) is T cxtrap 34, 
which is the same order of magnitude as time r ne q U ii ~ 60 for the system to reach metastable 
equilibrium, (b) Log- log plot of the equilibration time t s versus s for R = 20 and h = 1.258 and 
the same conditions as in Fig. |3(b)[ We find that r s ~ s x with the exponent x ~ 0.56. The 
extrapolated value of r s corresponding to the mean size of the nucleating droplet (~ 300 spins) 
is T ex trap ~ 30, which is comparable to the time r nequ ii ~ 40 for the system to reach metastable 
equilibrium. 

evolution of n^t), which can be fitted to the exponential form: 



We find that r s « 8.1 for s = 6. By doing similar fits for a range of s, we find that the time 
t s for the mean number of clusters of size s to become time independent increases linearly 
with s over the range of s that we can simulate (see Fig. EJ). The extrapolated value of t s 
corresponding to the mean size of the nucleating droplet (« 25 spins by direct simulation) 
is Tcxtrap ~ 34. That is, it takes a time of r extrap w 34 for the mean number of clusters whose 
size is the order of the nucleating droplets to become time independent. The time r ex trap is 
much longer than the relaxation time r s « 1.5 of the macroscopic quantities m(t) and e(t) 
and is comparable to the time r nequ ji ~ 60 for the nucleation rate to become independent of 
time. 

Because the number of clusters in the nucleating droplet is relatively small for R = 1 
except very close to coexistence (small h), we also consider a long-range Ising model with 



n s (t) = n S)00 [l - e~ t/Ta }. 



(4) 
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R = 20 and h = 1.258 (as in Fig. [3]). The relaxation time r s of the clusters near the 



pseudospinodal fits to a power law t s ~ s x with x ~ 0.56 (see Fig. 5(b)). We know of no 
theoretical explanation for the qualitatively different dependence f the relaxation time r s 
on s near coexistence (r s ~ s) and near the spinodal (r s ~ s 1 ^ 2 ). If we extrapolate r s to 
s = 300, the approximate size of the nucleating droplet, we find that the equilibration time 
for clusters of the size of the nucleating droplet is r cxtrap ~ 30, which is comparable to the 
time Tncquii ~ 40 for the nucleation rate to become independent of time. 

To determine if our results are affected by finite size effects, we compared the equilibra- 
tion time of the clusters for lattices with linear dimension L = 2000 and L = 5000. The 
equilibration times of the clusters were found to be unaffected. 

IV. STRUCTURE OF THE NUCLEATING DROPLET 

Because nucleation can occur both before and after the system is in metastable equilib- 
rium, we ask if there are any structural differences between the nucleating droplets formed 
in these two cases. To answer this question, we determine the nature of the nucleating 
droplets for the one- dimensional (ID) Ising model where we can make R (and hence the 
size of the nucleating droplets) large enough so that the structure of the nucleating droplets 
is well defined. In the following we take R = 2 12 = 4096, h = 1.265, and L = 2 18 . The 
relaxation time for m(t) is r g « 40, and the time for the distribution of nucleation times to 
reach equilibrium is r noquil « 90. 



We use the intervention method to identify nucleation [151 ] . To implement this method, we 
choose a time at which a nucleating droplet might exist and make many copies of the system. 
Each copy is restarted using a different random number seed. The idea is to determine if 
the largest cluster in each of the copies grows in approximately the same place at about 
the same time. If the percentage of copies that grow is greater than 50%, the nucleating 
droplet is already in the growth phase; if it is less than 50%, the time chosen is earlier than 
nucleation. We used a total of 20 trials to make this determination. 

Our procedure is to observe the system for a time £ b s after the intervention and determine 
if the size of the largest cluster exceeds the threshold size s* at approximately the same 
location. To ensure that the largest cluster at t Q b s is the same cluster as the original one, we 
require that the center of mass of the largest cluster be within a distance r* of the largest 
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cluster in the original configuration. If these conditions are satisfied, the nucleating droplet 
is said to grow. We choose t b s = 6, r* = 2R, and s* = 2000. (In comparison, the size of the 
nucleating droplet for the particular run that we will discuss is ~ 1080 spins.) 

There is some ambiguity in our identification of the nucleation time because the saddle 

n 

point parameter is large but finite [9J. This ambiguity manifests itself in the somewhat 
arbitrary choices of the parameters t Q bs, r*, and s*. We tried different values for t h s , r*, and 
s* and found that our results depend more strongly on the value of the parameter r* than on 
the values of t Q b s and s*. If we take r* = R/2, the nucleating droplets almost always occur 
one to two Monte Carlo steps per spin later than for r* = 2R. The reason is that the linear 
size of the nucleating droplet is typically 6 to 8R, and its center of mass might shift more 
than R/2 during the time t \, s . If such a shift occurs, a cluster that would be said to grow for 
r* = 2R would not be counted as such because it did not satisfy the center of mass criterion. 
This shift causes an overestimate of the time of the nucleating droplet. A reasonable choice 
of r* is 20% to 40% of the linear size of the nucleating droplet. The choice of parameters is 
particularly important here because the rate of growth of the transient nucleating droplets 
is slower than the growth rate of droplets formed after metastable equilibrium has been 
reached. Hence, we have to identify the nucleating droplet as carefully as possible. 

Because nucleation studies are computationally intensive, we used a novel algorithm for 



201 ] . The algorithm uses a 



simulating Ising models with a uniform long-range interaction 
hierarchical data structure to store the magnetization at many length scales, and can find 
the energy cost of flipping a spin in time 0((lnR) d ), rather than the usual time 0(R d ), 
where d is the spatial dimension. 

Figure [6] shows the fraction of copies for which the largest cluster grows as a function of 
the intervention time. For this particular run the nucleating droplet is found to occur at 
t « 37.4. 

We simulated 100 systems in which nucleation occurred before global quantities such 
as m(t) became independent of time, t < r g « 40, and 100 systems for which nucleation 
occurred after the nucleation rate became time independent (t > r ncq uii ~ 90). We found 
that the mean size of the nucleating droplet for t < r g is ~ 1200 with a standard deviation 
of a ~ 150 in comparison to the mean size of the nucleating droplet for t > r nequ ii of ~ 1270 
and a ~ 200. That is, the nucleating droplets formed before metastable equilibrium has 
been reached are somewhat smaller. 



10 



We introduce the cluster profile p c \ to characterize the shape of the largest cluster at the 
time of nucleation. For a particular bond realization a spin that is in the stable direction 
might or might not be a part the largest cluster due to the probabilistic nature of the bonds. 
For this reason bond averaging is implemented by placing 100 independent sets of bonds 
between spins with probability 

Pb = 1 - e- 2pj ^-^ in the stable direction. The clusters are 
identified for each set of bonds, and the probability pi that spin % is in the largest cluster is 
determined. The values of pi for the spins in a particular bin are then averaged using a bin 
width equal to R/4. This mean value of Pi is associated with p c \. Note that the spins that 
point in the unstable direction are omitted in this procedure. The mean cluster profile is 
found by translating the peak position of each droplet to the origin. 



Figure 7(a) shows the mean cluster profile formed after metastable equilibrium has been 



established (t > T ncqui i « 90). The position x is measured in units of R. For comparison we 
fit p c i to the form 

p(x) = Asech. 2 (x/w) + po, (5) 



with A c \ = 0.36, w c \ = 2.95 and po = by construction. In Fig. 7(b) we show a comparison 
of p c \ to the Gaussian form A g exp(— (x/w g ) 2 ) with A g = 0.35 and w g = 3.31. Note that 
Eq. gives a better fit than a Gaussian, which underestimates the peak at x = and the 
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FIG. 6: The fraction of copies for which the largest cluster grows for a particular run for a ID 
Ising model with R = 2 12 , h = 1.265, and L = 2 18 . The time for 50% growth is « 37.4. The 
largest cluster at this time corresponds to the nucleating droplet and has ~ 1080 spins. For this 
intervention 100 copies were considered; twenty copies were considered for all other runs. 
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(a) Comparison to Eq. ([5]). (b) Comparison to Gaussian. 

FIG. 7: Comparison of the mean cluster profile (•) in the ID Ising model after metastable equi- 
librium has been established with (a) the form in Eq. ([5]) and (b) a Gaussian. Note that Eq. (J3J) 
gives a better fit than the Gaussian, which underestimates the peak at x = and the wings. The 
x axis is measured in units of R. 



wings. Although Unger and Klein [12J derived Eq. (J5J) for the magnetization saddle point 
profile, we see that this form also provides a good description of the cluster profile. 

A comparison of the cluster profiles formed before and after metastable equilibrium is 
shown in Fig. [HI Although both profiles are consistent with the form in Eq. (j3J), the transient 
nucleating droplets are more compact, in agreement with the predictions in Ref. 1X9 . 



We also directly coarse grained the spins at the time of nucleation to obtain the density 



profile of the coarse-grained magnetization p m (x) (see Fig. 9(a)). The agreement between 



the simulation and analytical results 24J are impressive, especially considering that the 
analytical form is valid only in the limit R — ► oo. The same qualitative differences between 
the nucleating droplets that occur before and after metastable equilibrium is found (see 



Fig. 9(b)), although the magnetization density profile is much noisier than that based on 



the cluster analysis. 



V. LANGEVIN SIMULATIONS 



It is interesting to compare the results for the Ising model and the Langevin dynamics 
of the 4 model. One advantage of studying the Langevin dynamics of the <f) 4 theory is 
that it enables the efficient simulation of systems with a very large interaction range R. If 
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FIG. 8: The cluster profiles of the nucleating droplets formed before (dashed line) and after (solid 
line) metastable equilibrium has been established. Both profiles are consistent with the form given 
in Eq. ([5]), but the transient nucleating droplets are slightly more compact. The fitting parameters 
are A = 0.38 and w = 2.67 for the transient droplets and A = 0.35 and w = 2.95 for the droplets 
formed after the nucleation rate has become independent of time. 



all lengths are scaled by a large value of R, the effective magnitude of the noise decreases, 
making faster simulations possible. 

The coarse grained Hamiltonian analogous to the ID ferromagnetic Ising model with 
long-range interactions in an external field h can be expressed as 



(6) 



where (j)(x) is the coarse-grained magnetization. A dynamics consistent with this Hamilto- 
nian is given by, 



^ = ~M— + V = -M[-R 2 



2e(f) + 4w0 3 - h] 



(7) 



dt 5<t> dx 2 

where M is the mobility and rj(x,t) represents zero-mean Gaussian noise with 
(t](x, t)r)(x' : t')) = 2kTM5(x - x')5(t - t'). 

For nucleation near the spinodal the potential V = e<f) 2 + «0 4 — h<p has a metastable well 



only for e < 0. The magnitude of <fi and h at the spinodal are given by h s = a/ (8|£| 3 /27m) 
and 0s = (|e|/6it), and are found by setting V = V" = 0. The distance from the spinodal 
is characterized by the parameter Ah = \h B — h\. For Ah/h s 1, the bottom of the 
metastable well min is near </> s , specifically min = —<j) s (l + ^/2Ah/3h s ). 
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(a) Comparison with Eq. |[SJ). 
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FIG. 9: (a) The magnetization density profile of the nucleating droplets formed after metastable 



241 ] which has the form 



equilibrium has been established. The solid line is the analytical solution 
in Eq. (J5]) with the calculated values A = 0.085, w = 2.65, and po = —0.774. (b) Comparison of the 
density profile of nucleating droplets formed before (dashed line) and after (solid line) metastable 
equilibrium has been established by coarse graining the magnetization. The same qualitative 
differences between the nucleating droplets that occur before and after metastable equilibrium are 
observed as in Fig. [HJ although the magnetization density profile is much noisier than the cluster 
density profile. 



The stationary solutions of the dynamics are found by setting 5H/5(p = 0. Besides 
the two uniform solutions corresponding to the minima in V, there is a single nonuniform 
solution which approximates the nucleating droplet profile when the nucleation barrier is 
large. When Ah/h s <C 1, the profile of the nucleating droplet is described by Eq. (jHJ) with 

A = y/h s /6Ah/(J) s , w = (8h s Ah(j) 2 s /3)-^\ and p = min W- 



The dynamics ((7|) is numerically integrated using the scheme 25] 
<f>{t + At) = (f>(t) - AtM[ - R 2 —^ + 2e<p + 4u<p 3 - h] - 



At 



-'/• 



(8) 



dx 2 J ' V Ax 

where d 2 (j)/dx 2 is replaced by its central difference approximation. Numerical stability re- 
quires that At < (Ax/R) 2 , but it is often desirable to choose At even smaller for accuracy. 

As for the Ising simulations, we first prepare an equilibrated system with in the stable 
well corresponding to the direction of the external field h. At t = the external field is 
reversed so that the system relaxes to metastable equilibrium. We choose M = 1, T — 1, 
e = —1, u — 1, and Ah = 0.005. The scaled length of the system is chosen to be L/R = 300. 
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FIG. 10: Log-linear plot of the distribution P(t) of nucleation times for the one-dimensional 
Langevin equation with R = 2000 (x) and R = 2500 (•) averaged over 50,000 runs. The distribu- 
tion is not exponential for early times, indicating that the system is not in metastable equilibrium. 
Note that the nucleation rate is a rapidly decreasing function of R. 



We choose R to be large so that, on length scales of R, the metastable fluctuates near its 
equilibrium value min ~ —0.44. 

After nucleation occurs will rapidly grow toward the stable well. To determine the 
distribution of nucleation times, we assume that when the value of the field <p in any bin 
reaches 0, nucleation has occurred. This relatively crude criterion is sufficient for determin- 
ing the distribution of nucleation times if we assume that the time difference between the 
nucleation event and its later detection takes a consistent value between runs. 

Figure [10] compares the distribution of 50,000 nucleation times for systems with R = 2000 
and R = 2500 with Ax/R = 1 and At = 0.1. The distribution shows the same qualitative 
behavior as found in the Metropolis simulations of the Ising model (see Fig. [2]). For example, 
the distribution of nucleation times is not exponential for early times after the quench. As 
expected, the nucleation rate decreases as R increases. Smaller values of Ax and At give 
similar results for the distribution. 

To find the droplet profiles, we need to identify the time of nucleation more precisely. The 
intervention criterion, which was applied in Sec. IIVI is one possible method. In the Langevin 
context we can employ a simpler criterion: nucleation is considered to have occurred if <fi 
decays to the saddle-point profile (given by Eq. (JS} for Ah/h s <C 1) when is evolved using 
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FIG. 11: Comparison of the density profile <j)(x) of the nucleating droplets found by numerically 
solving the Langevin equation after metastable equilibrium has been reached for R = 2000 (x) 
and R = 4000 (•) to the theoretical prediction (solid line) from Eq. ([5]) using the calculated values 
A = 0.096, w = 3.58, and po = —0.44. The numerical solutions are averaged over 1000 profiles. 
The results suggest that as R increases, the observed nucleation profiles converge to the prediction 
of mean-field theory. 



noiseless dynamics 



261 ] . For fixed Ah these two criteria agree in the R — > oo limit, but 



can give different results for finite R 27]. 

In Fig. [TT]we plot the average of 1,000 density profiles of the nucleating droplets formed 
after metastable equilibrium has been established for R = 2000 and R = 4000. Note that 
there are noticeable deviations of the averaged profiles from the theoretical prediction in 
Eq. (jSJ), but the deviation is less for R = 4000. The deviation is due to the fact that the 
bottom of the free energy well in the metastable state is skewed; a similar deviation was also 
observed in the Ising model. We also note that the individual nucleating droplets look much 
different from their average. It is expected that as R increases, the profiles of the individual 
nucleating droplets will converge to the form given by Eq. ([5]). 

In Fig. [12] we compare the average of 1,000 density profiles of nucleating droplets before 
and after metastable equilibrium has been established. As for the Ising model, there are 



subtle differences consistent with the predictions of Ref. 



19l . The transient droplets have 



slightly lower background magnetization and compensate by being denser and more compact. 
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FIG. 12: The density profile of the nucleating droplets found from numerical solutions of the 
Langevin equation formed before (dotted line) and after (solid line) metastable equilibrium has 
been established. Nucleation events occurring before t = 15 are transient, and events occurring 
for t > 30 are metastable. Both plots are the result of 1000 averaged profiles with an interaction 
range R = 2000. 

VI. SUMMARY 

Although the time-independence of the mean values of macroscopic quantities such as 
the magnetization and the energy is often used as an indicator of metastable equilibrium, we 
find that the observed relaxation time of the clusters is much longer for sizes comparable to 
the nucleating droplet. This longer relaxation time explains the measured non-constant nu- 
cleation rate even when global quantities such as the magnetization appear to be stationary. 
By identifying the nucleating droplets in the one- dimensional long-range Ising model and the 
Langevin equation, we find structural differences between the nucleating droplets which oc- 
cur before and after metastable equilibrium has been reached. Our results suggest that using 
global quantities as indicators for metastable equilibrium may not be appropriate in general, 
and distinguishing between equilibrium and transient nucleation is important in studying 
the structure of nucleating droplets. Further studies of transient nucleation in continuous 
models of more realistic systems would be of interesting and practical importance. 

Finally, we note a subtle implication of our results. For a system to be truly in equilibrium 
would require that the mean number of clusters of all sizes be independent of time. The 
larger the cluster, the longer the time that would be required for the mean number to 
become time independent. Hence, the bigger the system, the longer the time that would 
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FIG. 13: The relaxation of the magnetization m(t) of the 2D Ising model at T = T c starting from 
T = 0. (a) R=1,T C = 2.269, L = 5000. (b) R = 128, T c = 4, L = 1024. The straight line is the 
fit to a power law with slope ~ 0.057 for R = 1 and slope ~ 0.51 for R = 128. 

be required for the system to reach equilibrium. Given that the system is never truly 
in metastable equilibrium so that the ideas of Gibbs, Langer, and others are never exactly 
applicable, when is the system close enough to equilibrium so that any possible simulation or 
experiment cannot detect the difference? We have found that the magnetization and energy 
are not sufficient indicators for nucleation and that the answer depends on the process being 
studied. For nucleation the equilibration of the number of clusters whose size is comparable 
to the size of the nucleating droplet is the relevant indicator. 

APPENDIX A: RELAXATION OF CLUSTERS AT THE CRITICAL TEMPER- 
ATURE 

Accurate determinations of the dynamical critical exponent z have been found from the 
relaxation of the magnetization and energy at the critical temperature. In the following 
we take a closer look at the relaxation of the Ising model by studying the approach to 
equilibrium of the distribution of clusters of various sizes. 

We consider the Ising model on a square lattice with L = 5000. The system is initially 
equilibrated at either zero temperature To = (all spins up) or at To = oo, and then 
instantaneously quenched to the critical temperature T c . The Metropolis algorithm is used. 

As a check on our results we first determine m(t) starting from T = 0. Scaling arguments 
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FIG. 14: The evolution of the number of clusters of size s = 100 at T = T c starting from To = 0. 
The fit to Eq, (H} gives n Si00 = 51.3, d = -42, C 2 = -15, n = 156, and r 2 = 1070. 



suggest that m(t) approaches its equilibrium value as [28] 



fit) = Bt~^ z + f 



(Al) 



where the static critical exponents are f3 — 1/8 and v — 1 for finite i? and /? = 1/2 and 
z/ = 1/2 in the mean-field limit. The fit of our results in Fig. [13] to Eq. ( IA1I) yields the 
estimate z ~ 2.19 for i? = 1 and z ~ 1.96 for i? = 128, which are consistent with previous 
results 29|, [30j. Note that no time scale is associated with the evolution of m(t). 

We next determined n s (t), the number of clusters of size s at time t after the temperature 
quench. Because all the spins are up at t = 0, the number of (down) clusters of size s begins 
at zero and increases to its (apparent) equilibrium value n Sj00 . The value of the latter 
depends on the size of the system. 

Figure [14] shows the evolution of clusters of size s = 100 for one run. Because we know 
of no argument for the time dependence of n s (t) — n Sj00 except in the mean- field limit 30], 



we have to rely on empirical fits. We find that the time- dependence of n s (t) can be fitted 
to the sum of two exponentials, 



n s {t) -n S)QO = C ie -' /ri + C 2 e~ t/T2 : 



where C\, t±, and r 2 are parameters to be fitted with r 2 > t±. 



(A2) 



Figure 15(a) shows the relaxation time r 2 as a function of s for R = 1 at T = T c starting 



from T = 0. Note that the bigger the cluster, the longer it takes to reach its equilibrium 
distribution. That is, small clusters form first, and larger clusters are formed by the merging 
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(a) T = 0. (b) T = oo. 

FIG. 15: The relaxation time r 2 versus the cluster size s at T = T c for R = 1 starting from (a) 
Tq = and (b) Tq = 00. The log-log plot in (a) yields r 2 ~ s 0A . 




(a) s = 30. (b) s = 3000. 

FIG. 16: The time dependence of the number of clusters of size s = 30 and s = 3000 at T = T c 
for R = 1 starting from To = 00. Note that n s= 3o monotonically decreases to its equilibrium value 
and n s= iooo overshoots its equilibrium value, (a) C\ = 2367, C2 = 332, n s= 3o,oo = 738, t\ = 16, 
and r 2 = 403. (b) d = -0.42, C 2 = 0.22, n s=300 o,oo = 0.11, n = 130, and r 2 = 1290. 

of smaller ones. The s-dependence of r 2 can be approximately fitted to a power law with 
the exponent 0.4. 

To prepare a configuration at T = 00, the system is randomized with approximately half 
of the spins up and half of the spins down. The temperature is instantaneously changed to 
T = T c . As before, we focus on the relaxation of down spin clusters. In contrast to the To = 
case, the evolution of the clusters falls into three classes (see Fig. [T6|) . For small clusters 
(1 < s < 40), n s monotonically decreases to its equilibrium value. This behavior occurs 
because the initial random configuration has an abundance of small clusters so that lowering 
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the temperature causes the small clusters to merge to form bigger ones. For intermediate 
size clusters (40 < s < 4000), n s first increases and then decreases to its equilibrium value. 
The initial growth is due to the rapid coalescence of smaller clusters to form intermediate 
ones. After there are enough intermediate clusters, they slowly coalesce to form bigger 
clusters, which causes the decrease. For clusters with s > 4000, n s slowly increases to its 
equilibrium value. The range of sizes for these different classes of behavior depends on the 
system size. In all three cases n s (t) can be fitted to the sum of two exponentials. One of 
the two coefficients is negative for 40 < s < 4000 for which n s (t) overshoots its equilibrium 



value. The relaxation time T2 is plotted in Fig. 15(b) as a function of s. 
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